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We present a fully non-perturbative determination of the 0{a) improvement coefficient csw in three-flavor 
dynamical QCD for the RG improved as well as the plaquette gauge actions, using the Schrodinger functional 
scheme. Results are compared with one-loop estimates at weak gauge coupling. 


1. Introduction 

Realistic simulation of QCD requires treating 
the light up, down and strange quarks dynami¬ 
cally. Incorporating a degenerate pair of up and 
down quarks have become almost standard by 
now, and a first attempt toward the continuum 
extrapolation has shown that the deviation of the 
quenched hadron mass spectrum from experiment 
[1] is sizably reduced [2]. Adding a dynamical 
strange quark is the next step, which has become 
possible with the recent algorithmic development 
for odd number of fermions[3]. 

The CP-PACS and JLQCD Collaborations 
have jointly started a 2-1-1 flavor dynamical QCD, 
employing the polynomial HMC (PHMC) algo¬ 
rithm for strange quark and the HMC algo¬ 
rithm for up and down quarks. We choose the 
renormalization-group (RG) improved action for 
the gauge fields, in order to avoid the lattice arti¬ 
fact present for the plaquette action [4]. We wish 
to use a fully 0(a)-improved action for quarks 
to control lattice spacing errors. Here we report 
on a non-perturbative determination of Csw for 
three-flavor QCD by the Schrodinger functional 
scheme both for the plaquette and RG-improved 
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Figure 1. AM as a function of M for the RG 
action with Nf = 3 at /3 = 2.2. 

gauge actions. 

2. Method and Simulations 

For the determination of cswj we basically fol¬ 
low the method of ref. [5], except for the choice 
B for the boundary weight of the RG-improved 
gauge action [6]. We refer to ref. [5] and refer¬ 
ences therein for notations in this report. 

We mainly use an 8^ x 16 lattice in our de¬ 
termination of Csw for tho RG-improved as well 
as the plaquette action with Nf = 3 dynamical 
quarks at several values of /3. Simulations with 
Nf = 4, 2,0 are also made for comparison. 
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We measure the modified PCAC quark masses, 
M and M', and their difference AM = M — M', 
at several values of csw and K. We have taken 
these parameters to realize M = 0 by an inter¬ 
polation, except at strong couplings for the case 
of -/V/ = 3, where an extrapolation to M = 0 is 
necessary as shown in Fig. 1. 

From the linear fit of AM as a function of M 
and Csw: = ao + aiM -l-a 2 CsW) we obtain the 

0{a) improvement coefficient csw = — 

ao)/a 2 , where AM^^'> = 0.000277, marked by the 
horizontal dotted line in Fig. 1, is the tree-level 
value of AM on a 8^ x 16 lattice. Note that the 
dependence of AM on csw becomes weaker at 
stronger couplings, so that the determination of 
Csw is more difficult, and hence the error is larger, 
at stronger couplings. 

3. Results 

In the upper plot of Fig. 2 we show the non- 
perturbative value of csw as a function of the 
bare gauge coupling for the RG-improved 
gauge action with Nf = 3(circles), 2(diamonds) 
and O(squares), together with the one-loop es- 
timate(solid line) and the mean-field(MF) esti- 
mate(dashed line) used in ref. [2]. Similarly, re¬ 
sults for the plaquette action with Nf = 3 (circles) 
and O(squares) are given in the lower plot of 
Fig. 2, together with the one-loop estimate(solid 
line) and the non-perturbative values by the Al¬ 
pha collaboration with Nf = 2(dotted lines) [5] 
and 0 (long-dashed line) [7]. 

In both cases, the non-perturbative values of 
Csw are almost Nf independent at weak coupling 
while they become larger for smaller Nf at strong 
coupling. This tendency can be clearly seen in 
Fig. 3, where csw is plotted as a function of Nf for 
the RG action(open symbols) and the plaquette 
action(solid circles). 

4. Comparison with perturbative esti¬ 
mates 

At first sight, the non-perturbative csw seems 
to undershoot the one-loop estimate at weak 
coupling for the RG action, while it converges 
smoothly from above for the plaquette action. 
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Figure 2. csw as a function of for the RG ac¬ 
tion (upper) and for the plaquette act ion (lower). 

We have found that the discrepancy seen for the 
RG action is caused by the one-loop 0{a/L) er¬ 
ror in Csw which becomes leading after the 
0{a/L) error at tree level is removed by requir¬ 
ing AM = AM°. In Fig. 4, the non-perturbative 
csw is compared with the one-loop estimate we 
have calculated on the same lattice size employed 
in the simulation, 8^ x 16. As seen from the fig¬ 
ure the non-perturbative value agrees with the 
one-loop estimate much better on the 8^ x 16 
lattice than in the infinite box. Note that the 
Oig^a/L) contribution to csw slightly depends on 
Nf through the fermion tadpole in the presence 
of the background gauge field of the Schrodinger 
functional scheme. Such an Nf dependence is in- 

^O(a) errors of cs^v in general cause 0{a?) errors in 
on-shell quantities, which are irrelevant in the Oia) 
improvement. 
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Figure 3. csw as a function of Nf for RG and 
P(plaquette) actions. 



Figure 4. csw for the RG action at weak cou¬ 
pling, together with the one-loop estimate on 
the 8^ X 16 lattice for Nf = 3(dotted line) and 
Nf = 0(dashed line). 

deed seen in the numerical data of Fig. 4. 

For the plaquette action we also confirm a pres¬ 
ence of the 0{gQa/L) correction, which is small 
on the 8^ x 16 lattice, as shown in Fig. 5 where 
csw is plotted as a function of o/L at (3 = 24 
with Nf = 0. The one-loop estimates (solid cir¬ 
cles) reproduce the non-monotonic behaviour of 
non-perturbative values(open circles) well. 

5. Discussion 

We have determined the non-perturbative 
value of csw for the RG action at several gauge 
couplings with Nf = 3,2,0. In order to obtain 
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Figure 5. Non-perturbative csw(open circles) 
and one-loop estimate(solid circles) as a function 
of a/L at /3 = 24 for the plaquette action with 
Nf = Q. 

an interpolation formula of csw as a function of 
(?Q, we have to eliminate large 0{gQa/L) correc¬ 
tions to csw present for the RG action. We are 
currently investigating this problem. 

We are also measuring the hadron spectrum for 
the RG action at /3 = 2 with Nf = 3 using the 
preliminary value of cswj in order to determine 
the corresponding lattice spacing. 
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